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ABSTRACT. We construct a nonstandard finite difference (NSFD) scheme for an SIRS mathematical 
model of respiratory virus transmission. This discretization is in full compliance with the NSFD 
methodology as formulated by R. E. Mickens. By use of an exact conservation law satisfied by 
the SIRS differential equations, we are able to determine the corresponding denominator function 
for the discrete first-order time derivatives. Our scheme is dynamically consistent with the SIRS 
differential equations since the conservation laws are preserved. Further, the scheme is shown to 
satisfy a positivity condition for its solutions for all values of the time step-size. 



1. Introduction 

In this paper, we will demonstrate how to correctly construct a finite difference scheme for a 
mathematical model of respiratory virus transmission [8 1 within the general discretization frame- 
work of Mickens' nonstandard finite difference (NSFD) methodology |E]|4j|5l. Our work clarifies 
the previous results of A. J. Arenas et al. JT|. While they began their calculations using a version 
of the NSFD discretization, they did not construct a scheme that is fully consistent with the NSFD 
methodology. In particular, they use an ad hoc (but excellent) selection for the denominator func- 
tions; further, their scheme places an upper bound on the value of the time step size that can be used 
for numerical computations. As we will show, our NSFD scheme has neither of these limitations. 

The following system ordinary differential equations have been used to model the transmission 
of the respiratory synacytical virus (RSV) transmission JTJ [8] 

(la) ^ = n - M 5(i) - p{t)S{t)I{t) + 7 J2(t), 

(lb) § = (3(t)S(t)I(t) - vl(t) - pl(t), 

at 

(lc) ^ = ul(t) - jR(t) - fxR(t) 

at 

where S(t), I(t), and R(t) represent the susceptible, infected and recovered populations, respec- 
tively. The initial conditions are 

(Id) 5(0) = So > 0, 7(0) = / > 0, R(0) = R > 0, 



2000 Mathematics Subject Classification. 39A50; 65L05; 65L12; 65L20. 

Key words and phrases. SIRS models; epidemic models; nonstandard finite difference schemes; nonstandard finite 
difference schemes; respiratory virus transmission. 

a The final version will be published in Journal of Difference Equations and Applications. 



1 



2 



RONALD E. MICKENS AND TALITHA M. WASHINGTON 



and the parameters (/j,, v, 7) are non-negative. The coefficient function f3(t) is taken to be periodic 
and non-negative, and often has the following general form UH) 

(2a) /3(t) =6 (l + 6icos(27rt + <j!>)), 

with 

(2b) < b , < 61 < 1, < <j> < 1, 

where is a normalized (constant) phase angle. Note that the structure of Equations ( flai l, ( flbb . 
and (TTcb is such that 

(3) {5(0) > 0, 1(0) > 0, P(0) > 0} {5(t) > 0, I(t) > 0, P(i) > 0}, 

for t > . We call this the "positivity condition." Adding Equations (flai l, d 1 bt . and (TTcb gives 
dP 

(4a) — = m(1 - P), P(0) = 5(0) + J(0) + P(0) = P 

where 

(4b) P(t) = S(t) + I(t) + R(t) 



is the total population. We call Equation (I4ab the "conservation law" corresponding to the system 
of the three ordinary differential equations given in Equation ([T). The solution to Equation d4ab is 


(5) P(t) = 1 + (P - l)e-"*. 

This result, along with the positivity requirement, allows us to conclude that if < Sq < 1, < 

I < 1, < R < 1, with < P < 1, then 

(6) < 5(i) < 1, < I{t) < 1, < R(t) < 1. 

Thus, we not only have positivity, but also boundedness of the solutions. (Note that P(t) is not 
necessarily a constant function which is what is stated in A. J. Arenas et al. iffl.) 

The differential equations (flai l. dTbl . (TTcb and (l4al > constitute a member of the class of population 
models satisfy a conservation law. Utilizing the procedures for constructing the appropriate NSFD 
scheme for such systems as shown in [5] gives 

t^-t N = hN, h = At; 

(S(t),I(t),R(t)) -> (S N ,I N ,R N ), 

P{t) /3(^) = fa. 

which is what is not implemented in the scheme shown in A. J. Arenas et al. fT). 

First, we observe that Equation (l4al > has the following exact finite difference scheme (HE] El 

(7) 7 = Ml - Pn+i), 



where the denominator function [2 1 is 

(8) </>(p, h) = 



efJ ,h 
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As explained in [5|, this denominator function must be used in the discretization of Equations ( TTab . 
( fTbl . and ( fTcb . otherwise, the exact NSFD scheme for the total population conservation law will 
not hold. In a similar manner, again see Mickens [5 1 for the details, the NSFD scheme for the three 
original differential equations, must take the form 



which is the exact finite scheme for the conservation law as expressed by Equation d4al >. Note that 
the solution for Pn is E| 



Thus, we can immediately reach the following conclusion: 

If, < S < 1,0 < I < 1, and < R < 1 with < P < 1, 

then < P N < 1, < S N < 1, < I N < 1, and < R N < 1 for all N > 1 . 

Consequently, this NSFD scheme gives numerical solutions that have the correct upper bound and 
they satisfying the positivity condition. 

The system of difference equations given in Equations (0 are not in a form suitable for com- 
putation. Since each of these equations is linear in Sn+i, In+i, and Rn+i, a rather long but 
elementary calculation gives J6) 





= - P N +i), Po = S + I + Rq, 



(10) 



P N = 1 + (P - l)e-" wv . 



(11) 




where 



(12a) 



A N 



e 3 ^ + <t>[{v + 7 ) + /3 N lN}e 2 ^ h + <t> 2 [iv + (v + y)p N I N ]e» h 



(12b) 



B N 



+ 



(lMf> + S N )[(1 + ij4) + u<f>][(l + n<j>) + 70] 
{"Iv^In + (70)[(1 + V<f>) + v4>)Rn 



(12c) 



Cn 



= {[(1 + IMj>) + 4>PnIn] [(1 + IMf>) + 7^] + {(flMRN 

+ (#jv)(m0 + S N )[(i + n4>) + j4>]}i n , 



(12d) 



d n = [(i + n4>) + <(>0 N i N ][(i + n4>) + "4>}Rn 

+ (0 2 ^ w )(^0 + S n )In + (WOK 1 + V<f>) + WnInVn, 
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and cj) = cf)(fi, h) is the denominator function as expressed in Equation ((8j. Inspection of Equations 
(fT2l shows that 

(13) A N > 0, B N > 0, C N > 0, and D N > 

for all N > 0, and for all step-sizes, h > 0. Thus, the above constructed NSFD scheme places no 
limitation on the step-size. 

We have achieved the stated purpose of this paper, namely, the correct construction of a NSFD 
scheme for an SIRS model of RSV. Our results demonstrate the power of the general NSFD 
methodology Q when applied to systems satisfying a condition of positivity and also having a 
conservation law Q. While this NSFD scheme may not be difficult to formulate, a successful 
implementation of a NSFD scheme requires that the construction of the scheme agrees with NSFD 
methodology. Our NSFD scheme has the following significant features: 

(1) The values of (Sjy+i, -fjv+ij -Rjv+i) are determined only by (Sjv, In, Rn), the step size 
h, and the non-negative parameters (/i, 7, v, bo, b\). 

(2) The denominator function is explicitly determined and does not require a particular spe- 
cialized form. 

(3) The scheme's solution satisfies the positivity requirement, the appropriate boundedness 
conditions, and maintains these properties for all h > 0. 

Finally, if we take as given the validity of the above stated SIRS mathematical model for RSV 
transmission HMD, then the NSFD discretization of this paper can be used to study in detail the 
dynamics of this epidemiological system. 
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